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ABSTRACT 


A general solution to the problem of one-dimensional thaw 
consolidation has been formulated by Morgenstern and Nixon (1971). 
In order to assess the validity of the theory it was necessary to 
develop a special oedometer (Permode) which could impose the 
necessary thermal and stress boundary conditions for one dimensional 


thaw consolidation. 


The permode permits the measurement of settlements, 
temperatures at various depths on the side of the sample, and excess 


pore pressures at the base of the sample during thaw-consolidation. 


Controlled thaw consolidation tests were carried out on 
three types of remoulded clays. The resulting data showed that the 
excess pore pressures and the degree of consolidation in a thawing 
soil depend primarily on the thaw consolidation ratio. The results 
obtained demonstrate that the theory adequately represents the soil 
behavior provided that the soil exhibits a linear relationship 


between void ratio and effective stress. 
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CHAPTER I 


INTRODUCTION © 
1.1 General Introduction and Review of Past Research 


In recent years there has been a greatly accelerated industrial 
development in the Arctic regions of Canada and Alaska. This 
development has been accompanied by an increase in the number of 
heavy structures being constructed on permafrost. The need for 
quantitative design criteria for foundations placed on permafrost 


has therefore become urgent. 


The basic information required in the design of any foundation 
is the strength and deformation properties of the subgrade soil. 
In general the strength and deformation of frozen soil is such as to 
present few foundation problems for the large percentage of structures 


where relatively light loads are involved. 


For foundations which thaw the soil as well as impose a load 
on it the situation can be much more critical. For structures such 
as large reservoirs, some knowledge of the pore pressures developed 
as thaw progresses is necessary in order to economically design 
the sand drains and embankment slopes. A knowledge of the rate of 
Settlement and the associated pore pressures is required in order to 
design foundatiuns for hot oil and gas lines constructed over 
perma-frost terrain. In general, the foundation design for such 
structures requires a determination of the following: 

1) The basic strength parameters for the thawed soil 
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2) The rate of are pressure development and 
dissipation as a function of the rate of thaw 
and the consolidation parameters for the soil, 

3) The deformations associated with the thaw- 


consolidation process. 


The basic strength parameters for a thawed soil can be 
determined by any of the test procedures presently employed in normal 


investigations of the mechanical properties of soils. 


The correlation between the rate of thaw and the associated 
pore pressures is of fundamental importance to the design of foundations 
and in assessing the stability of slopes. Tsytovich et al (1965) was 
the first to establish experimentally that the rate of settlement 
of a thawing soil was directly proportional to the rate of thaw. 
Once thawing was stopped, Tsytovich indicates that the rate of settlement 
would be aeeattien iy reduced. Tsytovich recognized two important 


aspects of the thaw consolidation process. They are: 


1) The pore pressures in the thawed zone dissipate 
sceundtine to the well known Terzaghi consolidation 
theory. 

2) The pore pressure slightly above the thaw boundary 
remains constant provided the rate of thaw is 


proportional to the square root of time. 


Tsytovich found that the pore pressures at the thaw boundary 


were some fraction of the applied load. That is: 
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U = (1 - m) BS (Oh) 


where U denotes the excess pore pressure at the thaw boundary, 
Be denotes the external applied load 


and m denotes a factor which is determined experimentally. 


Tsytovich provides values of m for different soils and 
formulates an equation for excess pore pressures using the above 
condition as a lower boundary. He was able to show that the excess 


pore pressures developed throughout the thawed portion of the soil 


depend not only on m but also on the ratio 


2VvC ; 
Vv 
where a denotes a tnermal constant indicating the rate of thaw, 
and — Cy denotes the coefficient of consolidation for the soil. 


It is difficult to assign a correct value to the coefficient 
m in design work and for this reason the pore pressures at the cee 
boundary were usually assumed to be equal to the applied load. This 
leads to over conservative designs in many cases. Furthermore, the 
coefficient gives us no indication as to what is physically going 


on at the thaw boundary. 


Zaretskii (1968) formulated a different boundary condition 
at the thaw boundary. He proposed that the change in water content 
of the soil as the thaw front moves a distance Ax is equal to the 


flow from the bcundary in time At. That is if: 
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where Aw denotes the change in water content, 
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S denotes the degree of saturation (= 1.0), 


Ae denotes the change in void ratio, 


and G. denotes the specific gravity of the soil grains. 
Then, 
6U 
Ve 
Nee SI 
G 5x (133) 
Yost 


The boundary condition at the thaw line can then be derived as: 


. G. cm ot (X,t) 
by fe U(X,t) = eee. 6x 
ot 


Zaretskii solves the classical consolidation equation subject 
to the above boundary condition. He is able to evaluate the coefficient 
m and therefore predict pore pressures and settlement. Zaretskii's 
assumption that the change in water content is equal to the volume 


of water which flows away from the thaw boundary is in error. 


W. G. Brown and G. H. Johnston (1971) attempted to show that 
the depth of thaw and the corresponding settlement can be predicted 
using transient heat conduction theory. The depth of thaw can be 
predicted fairly well in this way, however, there is no attempt to 
relate the rate of thaw to the consolidation parameters for the soil 
and therefore it is impossible to predict pore pressures or their rate 


of dissipation. 


A mathematical formulation for the thaw consolidation process 


has been derived by Morgenstern and Nixon (1971) using well known theories 
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of heat conduction and the classical Terzaghi consolidation theory. 
The approach is similar to that taken by Zaretskii but, there is a 
fundamental difference in the boundary condition employed at the thaw 


boundary. A summary of the theory is presented in Chapter II. 


1.2 Purpose of Thesis 


The purpose of this thesis was to ascertain the validity 
of the thaw consolidation theory as formulated by Morgenstern and 
Nixon (1971), by carrying out a series of temperature controlled one- 


dimensional thaw-consolidation tests on suitable remolded soil specimens. 


In order to carry out such experiments it was necessary to 
design a specially adapted oedometer (Permode), which could simulate the 
required temperature and stress boundary conditions. The initial 
temperature must be accurately controlled in the apparatus until some 
desired step temperature and total pressures are applied. Both heat flow 
and consolidation must be one-dimensional once the thaw-consolidation 


test begins. 


The permode was designed to meet the above requirements and 
permits the measurement of settlement, pore pressures at the base of 
the sample and temperatures on the side of the sample during the thaw- 


consolidation test. 


Comparison of the results of such tests with the predictions 
obtained from the Morgenstern and Nixon (1971) thaw consolidation theory, 
provide a good indication as to the validity of the theory. Such 
tests also serve to further the present understanding of the freeze-thaw 
phenomena of porous media and open up new areas for future 


investigation. 
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CHAPTER II 


THEORETICAL 


2.1 Summary of the Morgenstern-Nixon Thaw Consolidation Theory 


The mathematical formulation and solution to the problem of 
thaw consolidation which was derived by Morgenstern and Nixon (1971) 


is summarized here. 


A one-dimensional configuration is considered where a step 
increase in temperature is imposed at the surface of a semi-infinite 
homogeneous mass of frozen soil. Refer to Fig. (2.1). The solution 
to the heat conduction problem is attributed to Neumann and is given by 


Carslaw and Jaeger (1947). The movement of the thaw plane is given by 
X(t) = at!/4 (2.1) 


where X denotes the distance to the thaw plane 
and t denotes time 
a denotes a constant determined in the solution of the 


heat conduction problem. 


a = 2dvKo (ze) 
where Ko denotes the diffusivity of the thawed region 
and A denotes the root of a thermal equation which is a function 


of the initial temperature, the step temperature imposed, 
the conductivities of the thawed and frozen regions, the latent 
heat of fusion of the solid phase and the specific neat of the 


solid and liquid phases. 
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ONE - DIMENSIONAL THAW CONSOLIDATION 
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FIG, (2.1) 


In the thawed region, if it is assumed that the theory of 
consolidation is valid and that the total load does not vary with time 


then the governing equation is 


» dim. on 
V ae at (2.3) 
where Cy denotes the coefficient of consolidation 


x denotes the depth from the ground surface 


and u(x,t) denotes the excess pore pressure. 


The foregoing equations must be solved, subject to the following 
boundary conditions. 

(1) From equation (2.1) the initial condition is readily 
found to be: 

tae 0 eX = 20 

(2) Provided there is free drainage at the upper surface of 

the soil, 
one 6 SOR WS 05 aU 

(3) The boundary condition at the thaw line is derived as follows: 
The volume of pore fluid expelled from a small layer AX as the thaw line 
advances in an increment of time At is given by Darcy's Law, so that 

ou 
AV = - psa At (2.4) 
W 

where A denotes the cross-sectional area of a soil element 


k denotes the permeability of the soil 


SS (x,t) denotes the excess pore pressure gradient at the 


thaw interface 
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and tH denotes the unit weight of water. 


The discharge AV must be equal to the change in volume of a 


layer of thickness AX. The volumetric strain is therefore, 


ou 
av AV, K 3x (Gt) (2.5) 
VV AAX oe 
Yw dt 
But for a compressible soil, 
Nee 
View “mM Ao (200) 
where my denotes the coefficient of volume compressibility for the 
soil 
and Ac' denotes the change in effective stress. 
Substituting equation (2.6) into (2.5) will give, 
r SUD) 
ye Ae OX 
Ao' = Sem em AKe: (277) 
dt 
where Cos s : : 
Vesela ye the coefficient of consolidation. 
NE RLOCA SLY OSS saGaX eX eLS 
o(X,t) = aye + yX (2.8) 
where He denotes the stress applied to the surface 
and y devotes the bulk density of the soil. 
The pore pressure at x = X is, 
P(X, t) Seu Net) et Yy X (259) 
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The effective stress is thus 


o'(X,t)c= a ey eX wl X sit.) (2510) 
where y' denotes the submerged density of the soil. 
Now, 
Ao = ao" (X;t) - one (21,1) 
where a denotes the initial effective stress in the soil prior 


to any volume change on thawing. 


It is safe to assume that in ice rich soils on will be very 
nearly zero. Furthermore, if on is substantial because of stress 
nistory effects, it would contribute to stability and make less severe 
the problems that arise upon thawing frozen soils. From equations 


(2.10) and (2.11) with che = 0, we find 


Ao' = P, + y'X - u(X,t) (Oni) 


Substituting equation (2.12) into equation (2.7) gives the 
boundary condition at the thaw line as: 


et) 
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dt 
An analytical solution to the governing equations subject to 
the above boundary conditions has been found by Morgenstern and Nixon 


to be: 
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P 
u(X,t) =, erf( x ) * 1 
erf(R) + at eve yt ee 
vn R (2.14) 
where eas » and is called the thaw consolidation ratio 
oven 
Vv 
and , | 
erf( ) is the error function. 
ZV.CeC 
Vv 
Introducing the dimensionless variables 
eee : 1 
rz “tay ase the normalized depth of thaw 
and 'X(t) 
fee Tact , the ratio of the magnitudes of the applied load 
“ . 
to the overburden pressure, 
equation (2.14) becomes 
ie f= u(z,t) a 1 erf (RZ) . ee ees 
PRIy OX 1+W 2 1 
: Gert (R)is eg leg sit ye) 
2R r (2.15) 
vniR 


For a weightless material (W 


r 


= 0), equation (2.15) becomes: 


UZ Came erf(RZ) 
Re _pe 
erf(R) + = 
VaR (2.16) 


For a soil consolidating under its own weight (Wl, = «) 


equation (2.15) becomes: 
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The average degree of consolidation can be defined as the 
ratio of the consolidation settlement that has occurred at time t 
to the total consolidation settlement that would occur if thawing were 


stopped at time t. Thus, 


X 
a Peete Xo ez tL) tax 
S 0 
eas Uae 
$ =~ >= 
max 


X 
i le eye XU X 
fe) (e138) 


For the case of a weightless material (W, = 0) the average 


degree of settlement will be given by 


okey 
erf(R) + 
aero vr R 
ok’ 
(erf(R) + ) 
va R (2.19) 
and for vn = 0 
$=1-—y- 
lta 
2R (2-20) 


Values of excess pore pressures and the average degree of 
settlement for various thaw-consolidation ratios have been determined 


and are provided in chart form in the original publication. 
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2.2 Theoretical Predictions for the Oedometer Case 


It can be assumed that for the oedometer case Wh = 0. Thus 
the excess pore pressures will be given by equation (2.16). Values 
for excess pore pressure determined from this equation are presented 
in Fig.(2.2) The maximum theoretical pore pressures which will occur at 
the base of the sample in the oedometer for various R values can be 


derived from this figure and are presented in Fig.( 2.3). 


The average degree of consolidation for the oedometer case 
(W =e given by equation (2.19). The degree of consolidation for 


various R values is presented in chart form in Fig. (2.4). 


The pore pressures at the base of the sample in the post thaw 
phase of consolidation have been calculated by Nixon (unpublished) 
using a finite difference equation derived from the classical different- 
jal equation governing consolidation. The initial pore pressure 
distribution is taken as that which exists at the end of thaw and is 
therefore a function of the R value. The theoretical pore pressure 
dissipation curves evaluated for different thaw consolidation ratios 
is shown in Fig. (2.5). A theoretical post-thaw settlement curve can 
also be obtained for each R value and a set of such curves is 
shown in Fig. (2.6) This figure demonstrates that the theoretical 
settlement curve is almost insensitive to R for values of R less 
than 1.5. Furthermore it is apparent from this curve that if R is 
1.5 or less, the time factor at fifty percent consolidation has a 
value of 0.28 as compared with a value of 0.197 which applies for a 


uniform initial pressure distribution. 
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CHAPTER ITI 


DESIGN AND DESCRIPTION OF THE PERMODE 


3.1 Preliminary Testing 


a) Soil Description 

The soil used for the preliminary tests was Athabasca Clay. 
A summary of the mechanical properties of this soil are contained in 
Appendix A. The soil was selected because of its relatively high 
compressibility and low coefficient of consolidation (c = 15ft¢/year). 
Moreover, lean silty clays are more common in Northern Canada than 


clays with higher liquid limits. 


b) Preliminary Consolidation Tests 

Preliminary consolidation tests were carried out using standard 
consolidation apparatus with drainage from top and bottom. A slurry of 
Athabasca Clay was mixed with distilled water in an electric blender 
at a liquidity index of approximately 4.0. The slurry was warmed and 
deaired under a vacuum and stored in air tight containers for use as 


required. 


Three tests were carried out on samples which were initially 
consolidated under a load of 02kg/cm-. Load increments were added 
up to a maximum pressure of 6.3kg/em*. The ratio of 5 Was kept 


at a value close to 1.0 for all load increments. 


A second series of six tests was carried out using a slurry 


which was prepared as previously described. The slurry was placed 
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in cylindrical containers 2 inches long and 3 inches in diameter and 
was frozen in a deep freeze (-14°F). The mani were completely 
frozen in six to twelve hours. The samples were trimmed to size on a 
lathe while frozen and placed in the oedometer. The samples were 

then allowed to thaw at an uncontrolled rate in the oedometer under 
various initial loads. Once thawing and consolidation was complete 
under the initial load, subsequent load increments were applied up to 

a maximum of 6.3kg/cm2. In all of the above testing the load increments 
applied after the initial load were such that the load increment 


ratio was close to 1.0 


c) Test Results 

Values for tag for each load increment in the above testing 
were determined using the log time fitting method. A value for the 
coefficient of consolidation Cys coefficient of compressibility Mm.) 
and permeability k was computed from the consolidation data for each 


load increment. 


The values obtained for the three tests carried out on soils 
not subjected to freezing corresponded very closely to each other for 
each load increment. Typical values from one test are presented in 
Table 3.1. The corresponding values for a typical test carried out 


on the samples subjected to one freeze cycle are contained in Table 3.2 


The. e versus log o' curve for the tests is given in Fig(3.1). 
It is apparent from this figure that the application of a freeze thaw 
cycle to the soil tends to reduce the final void ratio obtained for 


the same load when the soil is subsequently thawed. That is, the soil 
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becomes over consolidated as a result of freezing. The difference 

_in void ratios is especially apparent at low pressures. The three frozen 
samples were thawed under different initial loads. It is apparent 

from Fig. (3.1) that the magnitude of the load on thawing has no 


Significant effect on the e versus o' curve obtained. 


It has been found that a plot of log permeability versus 
void ratio will produce a straight line for most soils. (Lambe 1951). 
Such a plot is presented in Fig. (3.2). Two straight lines were obtained, 
one for soils subjected to freezing and one for soils not subjected 
to a freeze cycle. It may be observed from this Frouteatnat for a 
given average void ratio, the permeability can be several orders of 
magnitude higher at low effective stresses, for a soil which has 


been subjected to freezing. 


This phenomena is evidently due to ice segregation during 
freezing which causes the formation of primary pore channels on 
thawing. As the slurry freezes, negative pressures are created at 
the penetrating frost line (Williams, 1967). Provided the rate of 
movement of the frost line is not too rapid, the ice will form nucleii 
and water will be drawn to the ice crystals. In this manner lenses of 
ice form leaving ribbons or pockets of soil in an over-consolidated 
state. The segregated ice was readily visible to the eye and was 
observed in all samples. When the sample is thawed there are, therefore, 
large primary pore channels left where the ice lenses had formed 
during freezing. It is these large pore channels which give the 
soil its very high permeability on thawing. I[t can be seen from 


Fig. (3.2) that as higher load increments are applied to the soil which 
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had been subjected to freezing, the primary pore channels are closed up, 
and the permeability approaches the values obtained for tests on samples 
not subjected to a freeze cycle. From Fig. (3.2) it can be seen that 
the permeability of the samples subjected to freezing is essentially 

the same as that found for samples not subjected to freezing at 

void ratio of about 0.75. This void ratio corresponds to an effective 
stress of approximately 1.5kg/cm°. This would indicate that the freeze 
cycle caused the soil to be overconsolidated to a pressure of about 


1.5 to 2.0kq/em?. 


Williams (1967) has shown that the magnitude of the negative 
pressures developed at a penetrating frost line depend primarily on | 
the air intrusion value of the soil, which in turn is a function of 
DaGtiC|ensi Ze. meuc mice | Ovdlt era tieverore,= tats tiesoverconsoy 1dau10n 
pressure achieved during freezing is a function of the air intrustion 


value. 


It is felt that further investigation should be carried 
out in this area to actually measure the magnitude of the negative 
pressures developed on freezing and relate these pressures to the 


overconsolidation pressure found on thaw. 
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3.2 PORE PRESSURE MEASUREMENT DURING CONSOLIDATION 


a) Design of the Pore Pressure Measuring System 


Before pressures could be accurately monitored in a thaw- 
consolidation test it was necessary to determine whether or not pore 
pressures could be accurately measured in a standard oedometer test. 
Numerous attempts to monitor pressures during an oedometer test 


have been made by others. The two major difficulties which must 


be overcome in measuring pore pressures accurately in the oedometer are: 


1) the effects of compliance of the measuring system 
on the pore pressures in the soil, and 


2) the effects of side friction. 


A fairly comprehensive theoretical study of the effect of 
compliance of the measuring system for the oedometer is presented by 


Perloff et al (1965). Perloff defines the parameter 


AR, 
rae) (351) 
where =A denotes the area of the sample, 


H denotes the height of the sample, 
my denotes the compressibility of the soil skeleton, 
A denotes the compressibility of the measuring system, 
and yn denotes the stiffness of the soil skeleton relative to 


that of the measuring system. 


He shows that the measuring system will have a minimal effect on the 


pore pressures being measured provided the ratio n is high. This is so 
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when the volumetric compressibility of the soil skeleton is high relative 
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to the volumetric compressibility of the measuring system. It is 
apparent from equation (3.1) that the effect of system compliance can 
be reduced by minimizing the volume of fluid in the measuring system 


and by increasing the volume of soil in the oedometer. 


Leonards and Girauld (1961) have investigated the effect of 
load increment ratio and side friction on pore pressures measured 
in the oedometer test. Their findings indicate that provided greased 
teflon is used on the sides of the oedometer and the load increment 
ratio is greater than a critical value (which varies depending on soil 
type and previous stress history), the pore pressures can be measured 
accurately and will conform fairly closely to that predicted by the 


Terzaghi consolidation theory. 


Barden and Berry (1965) investigated the effects of structural 
viscosity (creep) on the pressures measured in the oedometer tests. 
Their findings indicate that the effects of structural viscosity can 
be minimized by using thicker samples and applying load increment ratios 
greater than 1.0. Provided these conditions are satisfied, the pressures 
measured in a normally consolidated soil should be in good agreement 


with the Davis and Raymond (1965) non-linear consolidation theory. 


The permode (permafrost oedometer) was therefore designed to 
meet the requirements for accurate pore pressure measurement. A diagram 
of the apparatus is shown in Fig. (3.3). The pressure is measured 
through the use of a strain gauge which is connected to the transducer. 
Pressures can be measured to an accuracy of .002 psi. Originally a 
small porous stone was inset at the bottom centre of the permode. It 


was discovered, however, that better results could be obtained if a 
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circular porous stone which covered the entire base of the soil sample 
were used. The valves, poet ane’ as shown, permit the re-zeroing of 
the pressure transducer while a test is in progress, without affecting 
the pressures at the base of the soil. Re-zeroing was found necessary 
to compensate for changes in atmospheric pressure which occurred 


during longer consolidation tests. 


The load cap is fitted as loosely as possible without allowing 
extrusion of the soil around it. As shown in the diagram, a frictionless 
teflon guide is provided to ensure that the load cap travels smoothly 
downwards without tilting. This arrangement was found to give much 


better settlement and pore pressure data during testing. 
b) Test Procedure to Confirm Pore Pressure Measurement 


The permode was assembled with a rubber membrane stretched 
over the lucite ring which held it firmly in place. The inside and 
outside of the membrane were thoroughly greased. The measuring system 
was filled with deaired distilled water. No special procedures beyond 
the use of deaired water were found necessary to obtain good pressure 
response. A filter paper was placed over the bottom porous stone and 
the slurry was deposited in the apparatus under a head of approximately 
1/4 inch of water. The load cap was gently lowered into place and the 
sample was allowed to consolidate under the weight of the load cap 
(.02kq/em*). Once consolidation was complete the rubber membrane was 
released from the lucite ring. This procedure normally caused an_ 
increase in pore pressures which were allowed to dissipate before testing 
resumed. Various loads were then added to the sample up to a maximum 


2.16kg/em~. Time, deflection and pressure readings were taken as for a 
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standard consolidation test. Various load increment ratios were run 

on several different specimens. The load increment ratio was varied from 
0.8 to 11.4. Pore pressures were also measured after a specimen was 
unloaded, allowed to obtain equilibrium and then reloaded. The maximum 


overconsolidation ratio for this type of test was 4.0. 


c) Test Results 


It was found that the measuring system was sufficiently hard 
and that very nearly 100 percent pore pressure reaction could be obtained 
provided a sample of Athabasca Clay of approximately one inch in height 
was used, (volume 8lcc.). It is probable that for soils of lower 
compressibility a larger volume of soil would be required to obtain 


good response. 


Initially, oedometer tests were carried out without tne rubber 
membrane shown in Fig. (3.3). Typical pore pressure measurements 
obtained from such tests are shown in Fig. (3.4) and (3.5). The pore 
pressure curves were normalized to the theoretical Terzaghi curve by 


determining ts from the settlement data. Thus from 


fons 50 (203) 
ct 
VE Rey Ea 
and T=, = Sef t 
r 50 (3.4) 


It is apparent from Fig. (3.5) that at the higher loads there 


is a definite flattening of the pore pressure curves at the outset. 
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This flattening was found to disappear at the higher load increment ratios 


(oF S100). 


It was felt that the flattening may be due to friction and 
a greased rubber membrane was therefore fitted to the inside of the 
permode. Typical pore pressures measured with the rubber membrane in 


place are shown in Fig. (3.6), (3.7), and (3.8). 


The pore pressures shown in these figures conform very closely 
with the pore pressures predicted using the David and Raymond (1965) 
non-linear consolidation theory. It is of interest to note that when 
the rubber membrane was used without silicone grease, the flattened 


pore pressure curves as shown in Fig. (3.5) were again obtained. 


The pore pressures appear to be unaffected by the load increment 
ratio for the range of ratios considered, provided the rubber membrane 
is used. The pore pressure curves also agree with the David and 
Raymond theoretical pore pressure predictions for several tests performed 


on overconsolidated soils which had overconsolidation ratios up to 1.8. 


It was found that the single most important factor which 
affected the results obtained was side friction. [It is important that 
the load cap fit as loosely as possible, without allowing extrusion of 


the soil around it and the rubber membrane must be thoroughly greased. 


The results obtained confirm the accuracy of the pore pressure 
measuring system in the permode. It is apparent that further 
investigation could be carried out to determine the limiting conditions 
under which pore pressure measurements can be made accurately in the 


oedometer. Such factors as load increment ratio, overconsolidation 
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‘ia 
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P 
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Time Factor T = —— 
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Fig. (3.6) Base Pore Pressure Dissipation Curve for Standard Consolidation - With Membrane 


z ---- Theoretical Terzaghi 
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Fig. (3.7) Base Pore Pressure Dissipation Curve for Standard Consolidation with Membrane 
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ratio, volume of sample, soil compressibility and so forth will 


undoubtedly influence the pore pressures measured. 


3.3 Description of the Permode 


An overall schematic diagram of the thaw-consolidation test 
layout is shown in Figen eon The permode itself is shown in Fig. (3.3) 
and consists basically of a cylindrical lucite ring which is encased 
in a stainless steel outer ring. The lucite provides insulation and 
has a thermal conductivity of 5.0 x 103 cal/sec°C cc. which is an 
order of magnitude lower than that for most soils. It was anticipated 
that sampling of frozen soils could be most conveniently carried out 
with a sampler with a diameter of 3 inches. The inside diameter of 
the permode was therefore selected as 2.5 inches since this permits 
the testing of frozen undisturbed samples carved from a 3 inch diameter 
core. The length of the inside of the lucite ring is 3.5 inches from 
the porous stone to the top of the ring. This length permits the 
testing of high water content samples where large strains occur and 
also provides a more uniform thermal regime around the sample. The 
lucite is encased in a steel jacket 0.25 inches thick. This jacket was 
designed to limit radial deformation to 5.0 x 10 inches at an inside 
pressure of 20kg/em°. It was found that for most thaw consolidation 
testing, the pressure need not exceed 5kg/cmé and therefore, the 


steel jacket is not really required. 


The pore pressure measuring system has been described 
previously. Freezing of the distilled water in the measuring system 
during. thaw-consolidation tests is prevented by the addition of 50 
percent ethelene glycol. The lower porous stone is saturated with 


ethelene glycol before it is placed in the apparatus. 


The temperatures at each end of the sample are controlled 
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through the use of the thermo-electric elements located as shown, 

one at the base of the permode and one in the load cap. The temperature 
of these elements can be regulated within an accuracy of 0.2°C through 
the use of temperature controllers. The circuit diagram for a 
temperature controller is shown in Fig. (3.10). The power output of 

the controller is regulated by a thermistor sensor through a feedback 
system. The sensor is placed at the point where temperature is to 

be controlled and the specific temperature can be established by adjust- 
ing the potentiometer to the applicable resistance. It was found 
helpful to regulate the input voltage to the controller through the 

use of a Variac, in order to limit the current applied to the thermo- 
electric elements. The specifications for the thermoelectric elements, 
temperature sensors and the major components of the temperature 


controllers are contained in Appendix B. 


Four thermo-couples are located as shown, two on the hot 
sides of the thermoelectric elements and two on the top and bottom 
of the soil specimen, near the temperature controller sensors. The 
thermo-couples provide temperature readings to an accuracy of 0.2°C 
at two minute intervals on a Honeywell Strip Chart temperature recorder. 
Readings can be taken at various intervals down to a reading every 
fifteen seconds on this recorder. This system permits a continuous 
check on the temperatures at these points throughout the test. 
Additional temperatures are measured to an accuracy of 0.01°C through 
the use of thermistors located on the sides of the sample at two 
different heights as shown. These temperatures are recorded at one 


minute intervals on a digital-voltmeter recorder. 
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An adequate heat sink is provided for the lower thermoelectric 
element by placing it directly on a thin flat plate which is cooled by 
a refrigeration unit. During freezing of the sample a heat sink is 
provided for the top element by packing the top of the load cap with 
dry ice. It was found that the presence of air at room temperature 
around the top of the load cap during thawing of the sample provided 


a sufficient heat source for the top thermoelectric element. 


The sample was loaded directly through a counter-balanced 
load hanger. It was found that due to the large strains involved 
in testing frozen soils it was extremely difficult to maintain a 
constant load on the sample when levered loading systems were used 
since the lever arm had to be continuously balanced. A constant 
applied load is necessary in order to obtain accurate pore pressure 
data. The chief disadvantage of the direct loading system is that the 


magnitude of the applied loads is limited. 


Vertical displacements are read directly with the dial gauge 
or can be recorded on the digital recorder through the use of an LVDT. 
Both the dial gauge and LVDT are accurate to .0001 inch. A second 
dial gauge which is not shown in the figure was positioned so as to 
provide a direct reading of the sample height to an accuracy of 


.OO1 inches. 
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CHAPTER IV 
THAW CONSOLIDATION TESTS 


4.1 Test Procedure - Thaw Consolidation 


The Athabasca Clay was prepared as a slurry using the procedure 
described earlier. The permode was assembled with the rubber membrane 
stretched over the lucite ring. The pore pressure measuring system, 
the bottom porous stone, and the filter paper were saturated with 
ethyelene glycol. The temperature sensors were inserted to contact the 
membrane and the slurry was deposited in the permode under a head of 
1/4 inch of distilled water. The specimen was then consolidated under 
an initial load with the rubber membrane held stretched in place by 
the ring. Once initial consolidation was complete the rubber membrane 
was released and the extraneous pore pressures were allowed to dissipate. 


The specimen was now ready for testing. 


The reliability of the pore pressure measurements was 
normally corroborated by applying a nigher load increment and measuring 
deflection and pore pressure readings as for a standard oedometer test. 
Once pore pressures were completely dissipated, the transducer was cut 
off from the bottom porous stone with the valve and freezing of the 
sample commenced, using both top and bottom elements. Freezing of a 
one inch sample would normally take approximately two hours. Deflection 
readings were taken at approximately fifteen minute intervals. The 
sample was allowed to reach a constant temperature throughout (normally 


-10°C) and the next chosen load increment was placed on the hanger. 
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The top element was switched to heat and the step temperature selected 
(normally +10°C to +30°C). The step temperature would normally be 
achieved within three to six minutes, and would be accurately maintained 


by the temperature controller throughout the test. 


Deflection and temperature readings were taken at two minute 
intervals. This interval was found convenient and enabled the settlement 
and temperature curves to be accurately defined. Since the pressure 
transducer could be damaged by sustained negative pressures, it was 
necessary to check the pressure at intervals and leave the valve to the 
transducer turned off until positive pressures were obtained. Once 
positive pressures occurred, usually just prior to the thaw plane 
reaching the bottom of the sample, pressure readings could be 


recorded every two minutes. 


The time required to completely thaw a particular sample 
depends primarily on the sample height, the initial temperature, and the 
step temperature used in the test. An average length of time for thaw 
was one to two hours. Readings were continued after thawing was complete 
and the sample entered the post thawing consolidation phase. Once 
post thaw consolidation was complete and the excess pore pressures were 
completely dissipated, a standard consolidation test at a higher load 
could be carried out, or, the sample could be refrozen and another 
thaw consolidation test carried out. Once the testing program was 
completed for the sample, it was removed from the permode and the final 


height, moisture content and weight of dry soil were determined. 
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4.2 Calculations 


The height of the sample at any point in the test is found 
by adding the difference between the dial reading at that point and 
the dial reading at the end of the entire test to the height of the 
sample as measured when the sample was removed from the permode at the 
.end of testing. For the tests shown, the final height is the average 
of twelve random measurements on the sample. The void ratio at any 


times in the tests will be: 


7 H - ihe 
ate oa 
fe) et) 
wnere H denotes the height of the sample at that time, 
and Ho denotes the height of soil solids and is computed from, 
H apse 
Oo A Ge (4.2) 
where We denotes the weight of dry soil, 


A denotes the area of the oedometer (31.85 ena) 


> 


and G. denotes the specific gravity of soil solids (2.65). 


The degree of saturation of the sample at the end of testing 
can be computed from the moisture content determined at the end of 


testing. 


3 ase (4.3) 
where S denotes the degree of saturation, 
w denotes the moisture content of the sample at the end 
of testing, 


and e denotes the void ratio at the end of testing. 
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The various thaw consolidation parameters were calculated in 
the same way for each test. As an example, the calculations for test 
5 - 4 - 3 are presented. Refer to Fig. (4.1) and Table 4.1. A few 
words of explanation regarding the pressure, deflection and temperature 
data presented in chart form in Fig. (4.1) are necessary. The 
information has been plotted against the square root of time in all 
cases. The pressure is the positive pressure measured at the base 
of the sample expressed as a fraction of the total pressure on the 
sample. The height of the sample at the start of testing prior to 
freezing is shown and the point at which the post-thaw consolidation 
phase started is indicated by an arrow. 

H. denotes the height of the sample prior to freezing. 


He denotes the maximum sample height achieved when frozen. 


The void ratios at each stage of the test are given and 
are denoted as follows. 
e. the void ratio at the start of testing (prior to freezing). 
Cr the average void ratio of the sample once freezing 
js complete and the sample is at a uniform temperature. 
en the void ratio at the end of the complete thaw 
consolidation test. (Post thaw pore pressures 


completely dissipated). 


The temperatures were measured at the same points in 
all tests. e) 

T3 denotes the temperature at the top of the sample. 

Toy denotes the temperature at a point on the side of 


the sample a distance of 0.6 inches from the bottom 
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Ofethe sample. 

22 denotes the temperature at a point on the side of 
the sample a distance of 0.3 inches from the 
bottom of the sample. 


Ty denotes the temperature at the base of the sample. 
The following information is presented in summary form 
P. the pressure under which the sample was frozen. 


P. the pressure under wnich thaw-consolidation was 


carried out. 
P the load increment ratio = 


a the thermal constant for the test. 
c, the coefficient of consolidation as computed from 
the post thaw settlement curve. 


R the thaw consolidation ratio for the test. 


Sy The ratio of the settlement which occurred 


max during thawing to the total settlement which occurred. 


U., The maximum pore pressure obtained in the test 
pe 
0 as a fraction of the total pressure on the 
sample. 


From the settlement curve Fig. (4.1), the time at which the 
thaw phaSe reached the bottom of the sample is defined. It will be 


noted that the peak pressure was measured shortly after this time, and 
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that the temperature measured at the base of the sample was 32°F 

at this time. In this case, thawing was completed at /t = 10.45 min! /¢ 
or t = 109 minutes after the step temperature was imposed. This time 
becomes t = 0 for the post thaw consolidation process. The settlement 
for the post thaw phase can be plotted against the log of time on semi 
log paper. This nas been done in Fig. (4.2). The pore pressure 
dissipation curve is also shown on this figure. Once post thaw consolid- 
ation is completed the time for 50 percent post thaw consolidation is 


established as for a standard consolidation test. The construction 


is shown on the figure and value of teg = 13.5 minutes has been found. 


The coefficient of consolidation can be computed from 


As explained in Chapter II the te for the pressure 
distribution which exists at the start of post thaw consolidation has 
a value of 0.28 as opposed to a value of 0.197 for an initial uniform 


pressure distribution. 


The length of the drainage path is taken as the average height 


for the post thaw consolidation phase. From Table 4.1 H = 0.7/91 in, 


= 0:28 (.7791 x 2,54)% 
V 13 454% 60 
cv 1.355 x 107° Gina. 


An average value for m, can also be computed from the post 


thaw data. 
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Time - Minutes 
Post Thaw Settlement and Pore Pressure Dissipation Curve 


Test 5-4-3 


(4.2) 
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In the case of thaw-consolidation the increase in effective 
stress (Ac') must be taken as the total load applied to the sample. 
If Ao’ were taken as the load increment applied to the soil, the 
compressibility would have a value of infinity where no additional load 


was added , as for several of the tests. 
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The permeability k is thus 
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VeVeeW (4.6) 
290.28 (on) OO) we 
-/ 2 
Keene Sex| Omarcne/ sec 
The thermal constant a is determined from 
aX 
At es 

vt (457) 


where x is the average height of the sample for the thaw phase of the 
test, and t is the time at which the thaw phase reached the base of 


the sample. 


.8249 inches. Thus 


i 


From Table 4.1, x 


ce 0.8249 x 2.54 
10.45 /60 
a = .0257 anveccu& 


The value for R can now be computed. 
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Since the samples are frozen under a closed system, the degree 


of consolidation during thawing, Sy can be computed as follows: 


S 
max 
St = ies > 
nese SO Hac (4.9) 
where H. denotes the sample height at the start of testing, 


prior to freezing, 
Hot denotes the sample height at the time the thaw plane 
reaches the bottom of the sample, 
and Des denotes the sample height at 100 percent post thaw 


consolidation. 


Thus, in this example, 


>t = .8149 - 7947 
S—  {B149 =. 7635 
S 

oo = 0.39 

max 


A value for the maximum pore pressure U., is simply 
Po 


determined from the pore pressure data Fig.(4.1). 
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4.3 TEST RESULTS 


a) Test Series 5-4 


The soil used for this test series was Athabasca Clay. A 
summary of the mechanical properties of this soil is presented in 
Appendix A. A summary of the data obtained from this test series 
is contained in Appendix C, Table 1. The void ratio versus effective 
stress for the test series is shown in Appendix C, Fig. (C.1). The 


details for the stress path followed in this test series are as follows: 


Test 5-4-0 

A slurry at a liquidity index of 3.0 was placed in the 
oedometer and allowed to consolidate under a weight of .0894kg/cm®. 
The hanger counterbalance was adjusted to give an initial pressure 


reading of 100% of the total load on the hanger. Pore pressures 


dissipated completely within twenty-four hours. 


Test 5-4-1 

A load increment of .0845 was added to the sample and a 
Standard consolidation test was run to check out the pore pressure 
measuring system. A response of 98% was obtained and the pore pressures 
dissipated according to theoretical predictions. There was a significant 
time lag in achieving the initial 98% response. The settlement and | 


pore pressure curves are presented in Appendix C, Fig. (C.2). 


hest.5.—-4"="2 
The sample was frozen under the total load of 179kg/em? and 


a thaw consolidation test was carried out with no increase in pressure. 
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Unfortunately the lower thermoelectric element was left on and thawing 
did not occur proportional to the square root of time. The data for 


this test is presented in Appendix C, Fig. (C.3). 


Test 5 - 4 - 3 

The sample was frozen under the total pressure of 179kg/cm® 
and a thaw consolidation test was carried out under a pressure of 
0.367kg/cm?. The equipment performed satisfactorily. Data for this 


test is presented in Appendix C, Fig. (C.4). 


Test 5-4-4 
The sample was refrozen and a thaw consolidation test carried 


out with no increase in pressure. The data for this test is presented 


TBA OeNGIXe Cs er igne (Gs): 


rests 5. =)4-—- 5A 
A standard consolidation test was run on the sample with 
the pressure increased to 0.745kg/cm?. The data for this test is 


presented in Appendix C, Fig. (C.6). 


lest b= 4-5 


The sample was frozen under a pressure of 0.745kg/cm?. A 


thaw consolidation test was carried out under a pressure of 1.52kg/em?. 


Wheecata for this test 41s presented in Appendix €, Fig. (C€.7). 


Test 5 - 4 - 6 


The sample was frozen under a pressure of 1.52kg/em? and a 
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thaw consolidation test was carried out under a pressure of 2.220kg/cm*. 


The data for this test is presented in Appendix C, Fig. (C.8) 


When the sample was removed from the permode it was discovered 
that the tilt of the load cap was sufficient to allow a height difference 
of 0.15 inches across the sample. This situation was remedied by 
installing the teflon guide for the load cap as shown in the drawing 
of the permode. The time lag in achieving pore pressure response 
was also discovered to be due to the polyethelene porous disc used 
at the base of the sample which had a very low permeability. That 


porous disc was replaced with a more porous stainless steel disc. 
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b) Test Series 5 - 5 


The soil used for this test was Athabasca Clay. A summary 
of the data obtained from this test series is contained in Appendix D, 
Table 1. The void ratio versus effective stress diagram for this 
test series is presented in Appendix D, Fig. (D.1). The details for 


the stress path followed in this test are as follows: 


OS Byer Seer 

A slurry at a liquidity index of 3.0 was placed in the permode. 
The counterbalanced hanger was adjusted to give 100% pore pressure 
response under an initial load of 063kg/cm°. The pore pressures 


dissipated completely within twenty-four hours. 


WES BY Ty yz 

The applied load was doubled and a standard consolidation 
test was carried out. The pore pressures rose to 97% of the applied 
increment and subsequently dissipated very nearly as predicted 
by the Davis and Raymond theory. The settlement and pore pressure 


data for this test are presented in Appendix D, Fig. (D.2). 


est Se" 5 a9 3 
The slurry was frozen under the total pressure of 125kg/cm*.. 
A thaw consolidation test was carried out under a pressure of 0.250 


kg/cm*. The data for this test is presented in Appendix D, Fig. (D.3). 


During test 5 - 5 - 3 it was discovered that a faulty valve 


was allowing some anti-freeze to be sucked in at the base of the soil 
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during freezing. It is unlikely that this affected the results of 
test 5 - 5 - 3, however, the test series was terminated so that the 


valve could be replaced. 
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c) Test Series 5 - 6 


The soil used for this test series was Athabasca Clay. 
A summary of the data obtained from this test series is contained 
in Appendix E, Table 2. The void ratio versus effective stress 
for the test series is shown in Appendix E, Fig. (E.1). The details 


for the stress path followed in this test are as follows: 


Peste5 .-6.- | 

A slurry at a liquidity index of 3.0 was placed in the 
permode. The counterbalanced hanger was adjusted to give 100% 
pore pressure response under an initial load of .063kg/cm?. The pore 


pressures dissipated completely within twenty-four hours. 


Pest. 5-26 = 2 

The applied load was doubled and a standard consolidation test 
was carried out. PS wee pressures rose to 90% of the applied load 
increment. The pore pressures dissipated according to theoretical 
predictions, although there was some flattening of the curve which 
is believed to be due to friction. The settlement and pore pressure 


data is presented in Appendix E, Fig. (E.2) 


Test. 527-*67= 3 
The sample was frozen under the total pressure of 0.125kg/em®, 
and a thaw consolidation test was carried out under a load of 0.250kg/cm?. 


The data for this test is presented in Appendix E, Fig. (E.3). 
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Test 5 - 6 - 4 
The sample was refrozen and a thaw consolidation test was 


carried out under the same pressure (0.250kg/cm*). The data for this 


test is presented in Appendix E, Fig. (E.4) 
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The sample was refrozen under a pressure of 0. 250kg/cm” and 
a thaw consolidation test was carried out at an increased pressure 
of 0.935kg/cm°. The data for this test is presented in Appendix E, 
Bl Gea Eno) 


Test 5 - 6 - 6 


The sample was refrozen under a pressure of 0.935kg/cm", 


60 


and a thaw consolidation test was carried out at a pressure of 2.0kg/em®. 


The data for this test is presented in Appendix E, Fig. (E.6) 


fests o=—6"=— 7"=tO™o — 6 — 11 
The sample was refroZen and thaw consolidation tests were 
carried out under the same pressure (2.0kg/cm*) for five cycles. The 


data for each cycle is presented in Appendix E, Figs. (E.7) to (E.11). 


There is evidence from the temperature data dae Sine test 
series that the thermocouple at the base was not functioning correctly. 
The thermocouple read consistently warmer temperatures than were read 
on the other temperature sensors attached to the sample. It is thought 
that this may be the result of the thermocouple not being in proper 


contact with the soil at the base of the sample. 
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Hees te Se ese omen], 


The soil used for this test series was Athabasca Clay. 
A summary of the data obtained from this test series is contained 
in Appendix F, Table 2.- The void ratio versus effective stress for 
the test series is shown in Appendix F, Fig (F.1). The details for 


the stress path followed in this test series are as follows: 


ies Se dl 

A slurry was placed in the permode at a liquidity index 
of 3.79 (w = 95.7%). The sample was frozen under the weight of the 
load cap but freezing was carried out before significant consolidation 
occurred. Thaw consolidation was carried out under a load of .125 
kg/cm. A small amount of extrusion occurred during the initial 
stages of thaw-consolidation. The purpose of this test was to 
determine whether or not the theory would be valid for soils at very 
high water contents. The data for this test is contained in 
Ropengixers 1d. (bec). 
Test 5-7-2 

The sample was frozen under a pressure of 0.125kg/cm”, and 

2 


a thaw consolidation test was carried out under a pressure of 0.250kg/cm™. 


The data for this test is presented in Appendix F, Fig. (F.3). 


Ves Ge5u—3 (fe 3 
The sample was frozen under a pressure of 0.250kg/em?, and a 
thaw consolidation test was carried out with the pressure increased 


to 0.50kg/cm*. This test was an attempt to produce a very low R value 
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by having a very low step temperature increase on the sample (44.5°F). 
Under these conditions it was found that the time that the thaw plane 
reached the base of the sample was so poorly defined that the data 
could not be properly analysed. The data for this test is presented 


in Appendix F, (Fig. (F.4). 
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These tests were standard consolidation tests run on the 
sample. The settlement and pore pressure curves are not included 
in the data as the results were of little significance. The values 
for c.; k and my obtained from these tests are given in Appendix F, 


hab leat? 


| ofl nO nin. eteed nits nb ifnsno3’ deli! 
7 - 


Apa 
babrtant ton gm sitiho surgenth ayo hie 9 
2slifay: iF caret yhigta| staan! 30 ae mativedy! ale | 


7) i. } ak 
a x th ‘A ni avid ets etasd seni ay We 
s ope 


ite2 ee 
ven rae otal — on Ps 
eile. 4. 
AA eat bas a 


E A. 


a)» flestySeries! 5.-28 


The soil used for this test series was Lake Edmonton Clay. 
A summary of the mechanical properties for this soil is presented 
in Appendix A. A summary of the data obtained from this test series 
is contained in Appendix G, Table 1. The void ratio versus effective 
stress for the test series is shown in Appendix G, Fig.(G.1). The 


details for the stress path followed in this test series are as follows 


Test 5-8 - | 

A slurry was prepared, using the method previously described, 
at a liquidity index of 2.5 (w = 140%). The sample was frozen under 
the weight of the load cap, but freezing was carried out before 
Significant consolidation took place. Thaw consolidation was carried 
out under a load of .063kg/cm°. No extrusion was observed during 
thaw consolidation. This test was run to determine the validity of 
the Morgenstern and Nixon theory on a soil at very high water content. 


The data for this test is presented in Appendix G, (Fig. (G.2). 


Test to 8 <- 22 


The sample was frozen under a pressure of 063kg/cm*, and 


a thaw consolidation test was carried out under a pressure of 0.442kg/cm. 


The data for this test is presented in Appendix G, Fig. (G.3). 


The tests verified the theory for this soil, and it was felt 
pointless to continue with additional tests. The sample was therefore 


removed from the oedometer. 
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The soil used for this test series was Bentonite which 
consisted of 80 per cent montmorillonite and 20 per cent illite. 
All dissolved salts had nor washed from the soil and the sample 
had previously been consolidated to a pressure of 089kg/cm? 
and stored. A summary of the mechanical properties for this sojl is 
presented in Appendix A. A summary of the data obtained from this 
test series is contained in Appendix H, Table 1. The void ratio 
versus effective stress for the test series is shown in Appendix H, 
pide Hal). the details for the stress path followed in this test 


series are as follows: 


eSie5 295-5 

The sample was carved to fit the permode and was placed 
into the apparatus, being careful to eliminate trapped air. The 
pore pressures were permitted to come to equilibrium under a pressure 
of 063kg/cm~. The sample was frozen and a thaw consolidation test 
was carried out under a total pressure of 0.313kg/cm?. The data for 
this test is presented in Appendix H, Fig. (H.2). The post thaw 


data is presented in Fig. (H.3). 
Test 5-9-2 


The sample was again frozen under a pressure of 0.313kg/cm® 


and a thaw-consolidation test was carried out under a pressure of 
0.690kg/cm*. The data for this test is presented in Appendix H, 


Fig. (H.4). The post thaw data is presented in Fig. (H.5). 
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The tests on Bentonite were carried out primarily to 
verify the Morgenstern-Nixon theory at high R values. The results 
obtained verify the theory with the exception of the value for 

t obtained for test 5-9-1. The value obtained is considerably 


max 
higher than the value predicted by the theory. A moderate amount 


of extrusion was observed during the thaw phase of this test and 
is undoubtedly the cause of the high value of Pee obtained in 


this particular test. Max 
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CHAPTER V 


DISCUSSION AND CONCLUSIONS 


5.1 Validity of the Thaw Consolidation Theory 


Examination of the settlement curves for the thaw consolidation 
tests reveals that settlement is proportional to the square root of 
time provided the rate of thaw is proportional to the square root of 
time. This phenomenon was first observed by Tzytovich (1965) and is 
consistent with the theory developed by Morgenstern and Nixon. That 


is, the theory predicts that for a constant value of R 


< = const = J 

max (531) 
But, Swe ay Po X (572) 
So, Sy = J mM ian xX oe) 


The volume change associated with the ice-water transformation may be 


added directly to S 


t- 
So, So = J my) ve X + 0.09 1X (5.4) 
where So denotes the total settlement, 

and n denotes the porosity of the soil. 


From the above equation it is obvious that the settlement 
must be proportional to the depth to the thaw plane, which moves in 


this case at a rate proportional to the square root of time. 
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The Morgenstern and Nixon thaw consolidation theory also 
predicts that once thawing is stopped, the settlements will proceed 
at a rate predicted by the classical Terzaghi Theory using as initial 
values the excess pore water pressure distribution at completion of 
thawing. Examination of the settlement curves and pore pressure 
dissipation curves for the post thaw phase of any test reveals that 
this theoretical prediction is correct. 

U 

The maximum pore pressures obtained (5") have been plotted 
against the log R in Fig. (5.1). The agreement see the theoretical 
predictions and the experimental data is remarkable, however, the 


points are generally higher than the theoretical curve. 


It is also possible to plot the base pore pressure dissipation 
curves for each test. The real time is converted into time factor 


using the equation 


Saez st 
lea pes 
50 (5.5) 
wnere te is determined from the post thaw settlement curve. The 
values for 5 are then plotted against log T as has been done in 
fe) 


Fig. (5.2), for a sampling of R values selected from the test data. 

The majority of the dissipation curves obtained in testing were found 

to behave in much the same way as the curves shown in Fig. (5.2). It 
will be observed from this figure that although the shape of the 
experimental curves approximates the shape of the theoretical curves, 
the experimental curves are shifted to the right. This same phenomenon 
is observed in pore pressure curves obtained from standard consolidation 
tests, where the actual pore pressure curves are shifted to the right 


of the theoretical Terzaghi curve. It is possible that the shift of 
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the dissipation curves found in post thaw consolidation may be due 
to nonlinearity of the void ratio versus effective stress 


relationships. 


The degree of consolidation a has been plotted against 
log R in Fig. (5.3). The agreement between the experimental data 
and the theoretical prediction is only fair in this case. Possible 
reasons for the scatter of points on this figure will be discussed 


later. 


Examination of the e versus log o' curves for test series 
5-4 and 5-6 reveal that the void ratio is significantly reduced as 
successive freeze-thaw cycles are applied even though the load on the 
Sample remains constant. This is probably due to local overconsolidation 
which takes place during the freezing process when high negative 
pressures are generated in the soil. Presumably the sample swells on 
thaw but does not return to its original void ratio overall. There is, 
therefore, excess water available, which flows out of the sample 
according to the standard consolidation process. The phenomenon explains 
why the pore pressures on thawing can be higher than the load increment 


imposed an the test: 


It is apparent, therefore, that the value of on wcune 
effective stress prior to consolidation, must be a function of the 
pressure to which the soil is overconsolidated during freezing and the 
coefficient of evenin for the soil. In the pore pressure calculation 


for these tests it was assumed that oh was zero in all cases and hence 


all pore pressures were expressed as a percentage of the total load on 
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the sample. The data shows that this assumption was essentially correct 
for the soils tested. It is conceivable, however, that for soils 
subjected to freeze-thaw cycles where the negative pressures generated 
during freezing are insufficient to cause a high degree of overconsol- 
idation, or where the coefficient of swelling is very nearly equal 

to the coefficient of consolidation, then Te could have a magnitude 
almost as high as the final effective stress and hence the pore pressures 


achieved on thawing would be much lower. 


The prediction of pore pressures during thaw consolidation 
requires, therefore, a better understanding of the negative pressures 
generated during freezing and the role these pressures play in 
determining the initial effective stress on thaw. Considerable work 


has been carried out by Williams (1967) and others in this area. 


Special attention should be paid to the data obtained from 
Pestse oe /P el Appendix ts nig-.(fec)sendetest 5-6 — | Appendix’ G, 
Fig. (G.1), which were thaw consolidation tests carried out on 
soils which were initially at very high void ratios. It will be 
noted that despite the large strains involved, (approximately 50% 
strain), the resulting settlement curves display the general form 
predicted by theory. The fact that this should occur is significant in 
that the theory assumes overall strains are small. The values for ae 
obtained from these tests agreeSreasonably well with theoretically Hers 
predicted valuesUFig. (5.3). From Appendix F, Fig. (F.2) it will be 
seen that the pore pressures behaved quite erratically just after the 
maximum pore pressure was read. The phenomena was observed to occur 


to a very minor degree in several of the other tests run at lower 


void ratios. It is thought to be due to the collapse of the soil 
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Structure in a localized pocket where a large ice crystal had formed. 
Such a collapse would naturally be accompanied by a rapid settlement 


and a significant rise in pore pressure. 


The maximum pore pressure achieved in test 5 - 7 - 1 was 
99.2% of the total load. This pressure is much higher than the 
theoretically predicted value of 29% based on the R value of 0.431. 


There are several possible explanations for this discrepency. 


FAVS Ofealle, theavalierior Cy which was used to calculate 
R was determined from an erratically behaved post-thaw settlement curve 
Appendix F,°Fig.s (F.2). This results in a value for Cy which is not 
completely reliable. The values for C,, have been plotted against 
void ratio for all the thaw consolidation tests performed on Athabasca 
clay (See Fig. 5.4). The variation of c, with void ratio has also 
been plotted for two other soils, Chicago clay (Taylor, 1948) and 
Tyee River clay (Chattopadhyay, 1969) and these curves provide a 
better indication of the behaviour one can expect in C, at different 
void ratios. The points obtained from thaw consolidation tests on 
Athabasca clay show considerable scatter, however, a fairly definite 
trend results as indicated by the curve. The trend indicates that Cy 
should Yemain very nearly constant as void ratio increases, however, 
further tests should be carried out at high void ratios to confirm this. 
The value for c, obtained from test 5 - 7 - 1 (circled) is clearly 
anomalous and sh uld have a value of approximately 1.0 x 10°? cm/sec. 
Pamuli smc ccCi mye UCmiomUScdaLoucal culate Rethensayvalue Tor ReoT 1935 
is obtained. This value for the thaw consolidation ratio predicts 


a maximum pore pressure at the thaw plane of 75% of the total load 
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Fig. (5.5) Normalized Depth of Thaw Plane versus Normalized Square Root of Time 
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which is closer to the experimental value obtained. 


The high pore pressures obtained for test 5 - 7 - 1 might 
also be due to the non-linearity of the e - o' relationship which is 
known to be increasingly non-linear at lower effective stresses. For 
very high ice content soils, in fact, it is probable that considerable 
change in void ratio would have to take place before there was any 
increase in effective stress. Such a condition would act to maintain 
pore pressures for longer periods than predicted using a linear e - o' 
model. It is important to note that the present theory is in error 


on the unsafe side in this respect. 


It is strongly recommended that several thaw consolidation 
tests be run on-soils at very high void ratios before any definite 
conclusions are drawn regarding the validity of the Morgenstern-Nixon 


thaw consolidation theory for soils under these conditions. 
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Dace. SourcesroTeerror 


The relative deflection can be measured sufficiently accurately 


with the dial gauge where overall settlements (S._.) are at least 


max 
0.01 inches. For tests where the overall settelement is less than 


this value, the dial gauge positioned as it is to read deflections 


of the load hanger, did not give reliable readings. It is felt that 
S | } 

the scatter of results on the < versus log R curve is due partly 
max 

to the unreliability of the dial gauge readings. This is especially 


true where the thermal component of settlement becomes large with 
respect to consolidation settlement. For example, in test 5 - 6 - Ill, 


Sy was found to be 0.0007 inches and Sek 0.0045 inches. A very small 


vibration or movement will change the ratio = by a substantial 
max 


factor. 

It was also discovered that temperature changes above the 
freezing point in the soil will cause the height of the sample to 
change. As the temperature dropped from room temperature (75°F) to near 
freezing, the height of the sample would decrease by as much as 0.01 
inches. This decrease in height was always accompanied by a small 
negative pore pressure which was measured at the base. As the pore 
fluid began to freeze the negative pressures increased and of course, 
the sample height again increased. The temperature effects on the sample 
height will affect the value of *t to a greater or lesser extent 


max 
depending on the overall strains involved in the particular thaw 


consolidation test. 


The pressures measured in the thaw consolidation test are 
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quite accurate. The process of verifying the accuracy of the pressure 
measuring system at the start of each test series by carrying out 

a standard consolidation test ensured that the system was functioning 
correctly. It was found that the pressures measured were affected 
slightly by the temperature of the soil. As the temperature of the 
soil rose in the post thaw phase, the pressures measured would tend 

to rise. The rise in pressure counteracted the pressure drop due 

to dissipation to a small degree. The effects of temperature on 

the excess pore pressures obtained were never more than 1 to 2% 


of the total stress on the sample. 


Every attempt was made to freeze the soil as quickly as 
possible and to do so under closed system conditions. The sample 
was frozen from the top and generally froze completely within thirty 
minutes after the freezing plane entered the top of the sample. 
There was undoubtedly some migration of water from the bottom of the 


sample to the top during the freezing process but the amount would 


be negligible due to the low permeability of the soil (1 x 1084 cm/sec). 


By closing the lower porous stone to the atmosphere and freezing the 
sample at the top with dry ice, the soil could be completely frozen 
with virtually no change in overall moisture content. The void ratios 


which are shown in Table 1 of Appendices C to H for the frozen soils 


are, in most cases, 10% higher than the void ratio which existed before 


freeze commenced. This confirms that there was no substantial increase 
in overall water content as a result of the freezing process and that 
the increase in void ratio was due almost entirely to the phase change 


of water to ice. 
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The temperatures and step temperature increases were 
accurately maintained by the temperature controllers and no problems 
were encountered with this system. The recorded temperatures are 
accurate, but it proved difficult to determine precisely where the 
thermocouple or thermistor is measuring temperature. The thermistors 
on the side of the sample could easily to 0.05 inches higher or lower 
on the side of the sample. The thermocouple at the base measured 
consistently higher temperatures in test 5 - 6, than were expected. 


This thermocouple is thought to be defective. 


Although no direct observations could be made as to the 
progress and shape of the thaw plane, it is possible to examine the 
available data to determine whether or not the thaw plane was planar 


and moved downwards at a rate proportional to the square root of time. 


Examination of the settlement curves for tests 5 - 6 - 8 
to 5 - 6 - 11 where 90% of the settlement was due to the phase change 
of ice to water, reveals that during thawing, settlement occurred 
proportional to the square root of time. Assuming the porosity of 
the soil is constant it is apparent from these curves that the thaw 


plane must have moved at a rate proportional to the square root of time. 


The progress of the 32°F isotherin as it moves downward 
through the sample can be plotted from the temperature data. The data 
for test series 5 - 4 and 5 - 7 has been plotted in Fig. (5.5). .Data for 
test series 5 - 6 is excluded because of the unreliable base temperature 
readings. The height of each thermistor has been normalized to the 
average height of the sample for the thaw pnase of each test. In this 


case, since the temperature at the top of the sample was measured at 


- a 
~ om fda 


oni? 


63h 


907 fi) 


# 


a ; 
eentw auc er dl sae wre ete oF 
as 
nq nn ber PelliMymnos ai 


eorsorhiied Mebane BAT Si tive 


Mab Zac 5 GE igang! aaa eta Anite ial | 
“ae orintw pagent ary ree banat vain 0066 
$i) 07H nh htSoty) nant eo bs Pag 
ow Posnges ilies | 
Sreit 
te uy. seh iy 


_ 1 ; 


— 7 


_ 
4 3 
ww 
e 
» Ah 
i = 


ororntirsteh of oar 


att ‘Jowinhee(aimsd pa Pivecian ot 1 nat rie 
* shitoat 20.0.6) Uf Meteaubiues on 
> . Pa q ‘ . 


a 

. =) Ayu 

ow Hen ott Se Anqyvooon{Sns Si Tar oat Se of 
, . n il : 

. ‘ x 5 7 7 

cat ) - desl Wr Vout evegiigy istp  Yaines2 feng; 


Tost yb id’ ir Sy y i ef tay 2 ott. 2: AT 
2 FtogT yh ag eT Sy, is tq0 cH gxi rf 


9 . : _ 7% | . 2 ii = oars 
gh eK ad Li luod Zhi BvIReIh Coane On Hove TA =; 
2: Se ie 

t coneiqeweay, ody ta Sante 5 chutes 
i a i \ ai tie rom ye J 3 i Th igsoh 8 sey of és io! ie) 
ed? of lenohivegoig afi 69% sinned 


are. OP ee 
us: ‘ ea Aoniylhifog-ans +o. Hor tonne’ 
A] as 


t invstis:! a Ad pain ons 2laavers ener nf Yo 


Thi 43! , fa 10 “ | ree 4 iJ Ce tina aie “ rq 


f bh f: Aart ta > oy. 3 ) UY Sn: 25 -@ o3 


1 B28). WOT) SEARS 2b 0g *nede ae ah a 

s 7 r 

pipe 2 6) [enohh(oqaid Sher sore Ge aon il eum ace Te 
> ‘ ’ ' a a 


- _ 7 ’ _ 
é 


aid ead Moat i. 
acd 7" _ ae 

tis ‘Seatig He oi din ali 
eae 


shat. cue by eusyt st o eer 


a 
sn 3 - i Hod aft ‘ou iyi 


a 7 


fe) 


the top of the upper porous disc, the average height of the sample 
includes the thickness of the porous disc (.125 inches). It is assumed 
that the thaw plane moves through the porous disc at the same rate it 
moves through the soil. The value for the square root of time at which 
the 32° isotherm reached each thermistor has been normalized to the 
value for the square root of time when the 32° isotherm reached the 
base of the sample for each test. A plot of normalized depth versus 
normalized square root of time will, therefore, be independent of the 

a value for each test and the points should lie.on the dotted line 
shown in Fig. (5.5) provided the 32° isotherm moves downwards 
proportional to the square root of time. Since the thermistors are 
located on opposite sides of the sample the points will be on the 
theoretical line only if the 32° isotnerm moves downwards in a plane 
which is very nearly horizontal. The data plotted on Fig. (5-9) 
confirms that the 32° isotherm did move downwards in the permode 
proportional to the square root of time. The scatter of the points 

is probably due to the fact that the height of the thermistors could 


not be precisely determined. 


The base of the permode should simulate the existance of 
an infinite depth of frozen soil below the sample itself. In order 
to check this condition, a computer program was run to determine the 
temperature isochrones for a semi-infinite soil mass subjected to the 
step increase in temperature which was used in the test. The solution - 
employes a finite difference method formulated by Ho et al (1970) and 
was implemented on the computer by Mr. J. F. Nixon. The values for a 
and the moisture content of the soil as well as a typical unfrozen 


moisture content versus teinperature curve for the soil were incorporated 


ut barvoys” a3 abi 


7 ~_ Ms Sa 
Petes 2 hE (ositink: ite). 23) i nat — 
: Co We 7 — 7 _ eee 


- 


. ; a ; 
| DY Se emia oft te were prey att Futter 
dxidw 26 ate Jo dom & sae gilt “ot autew BAT... 
= ; ‘2 & : 7 
eat. ah beni lotro ansed 2a vers liisty Agee: 41 MIUNIG 
ie tr baiscety 0r Parag OOS ely rey Wi S004 mee 3d i 
ene wig), bowl scrron Wd Jools A .deos- Nee 4or. ofqmna ott Ye — 
ari i? Ayin? fi hctit Sal a Tsk} 775 Aj i rie ar.) ot $004" siya 
sei! beddon ofl) | Sivas Fog. uth. bite) Feed eg an. svt 
: chiara eovan mBAsoa OLS od, Hab Pred f0bd) oo mt wore 
; : i - _ 
V6. 2iolejmie ans 2Hre “92.70 Soot eTilina ont Ge nest svoqing, 
ais no ol TTI mrad ery « 1 Wy te pais af iia 0 Wo bs s507 
. | . ; | . : ’ a 
"M6IGQ B WT Veen ‘on wyandoer’° Nee art Th Ue ont] thates ost: 


(F-2) et? no tiniacty agph Sit” Témnestion Y ta ci of wolite 


ts P 
ebourray | | CH bei av euiseph OSE aged iy ‘win tnoD 
“ ‘ ? , 
2 P ‘ 4 feak o4 = 
¢-i1 vy iy ‘ 4 oh. a0C4. ey r ne. ot isis 
‘a ; ; = aa a al 
te pee mn 7 , 
Big iar Ly 


7 iow? Sig 6? aub cideamg at 
: . : > = 

ven linetah \ far toang og Jon 
: ry 


a 
“ : aa 
’ 4 : Py : 4 i i : - 
F 1 =e hit ev AS : : li Peat - - in ae 2a sit _ 
m ns 2, 


wre Wt «Teas 7 siqpor ely woled! (nz (eset Wages = nee 


eft) ofumeset ot pon 2 va Ty ao ig veane Ih mola ald + 


fais: of 453 2gkdus 2 | wih Eloe sia ve oe 1 eaniortg a ie ag 

nok wiwe ath i aa OF; 5rid nm b she Gi" sts cons f bs 
; ic 

baw (OTR) De oh ‘baaslame 


ei is 


ieee Pe ta. a 
—_ _ 
n 


te 


ee hs 
Mesracc’ .atay) agef fs v" sat touz 4 


aig iii Sha , Haowit om rian 
4 } 1 : : 


80 


into the solution. From the temperature isochrones derived in this 
manner it was possible to determine what the temperature at different 
times should be at a depth equal to the height of the soil sample. 

A theoretical base temperature versus time curve can thus be established. 
Such a curve is shown in Fig. (5.6). The actual points as measured 

by the thermocouple at the base of the soil for Test 5 - 4 - 6 have been 
plotted on this figure. The agreement is very close and is further 
evidence that the assumed heat boundary conditions were being 

correctly simulated in the permode. The base temperatures for all 

tests were checked in this way and very close agreement was found except 
fOracest series 5 — 6. Again, this is, in all probability, due to 


the faulty thermocouple used at the base for this test series. 


The temperature data can be plotted in the form of temperature 
isochrones for any particular test. As an example, the temperature 
isochrones for test 5 - 4 - 4 have been plotted in Fig. (5.7). The 
Shape of a typical isochrone determined from the computer solution 
previously described, has been plotted on the figure for comparison. 
The isochrones derived from the experimental results demonstrate a 
Shape which is close to that predicted by the computer solution. The 
theoretical solution predicts a sharp break in the isochrone at the 
32° isotherm as shown on the diagram. It is apparent from the 
experimental data that a sharp break does not occur, rather, that there 
is a smooth curve at this point. It is likely that this deviation 
is due to the fact that the thaw plane is in reality a zone of thaw 


rather than a planar surface. 


The time that the thaw plane reaches the bottom of the sample 
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Temperature at Base of ‘Soil °F 
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Theoretical Base Temperature 
2 Measured Base Temperature - 


(Tie - min! /? 


. (5.6) Variation of Temperature at Base of Sample with 
Square Root of Time 
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Fig. (5.7) Temperature Isochrones for Test 5-4-4 
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can be either well defined by the temperature, deflection and pressure 
data as for test 5 - 4 - 3, or very poorly defined as for test 5 - 4 - 6. 


The value for R is mildly sensitive to the selection of the time of 
5 

end of thaw, however, the value for < is very sensitive to the height 
max 

corresponding to the time of end of thaw. Some attempt was made to 


standardize the time selected as the end thaw time. In general the 
dial reading just prior to the maximum pore pressure reading seemed 
to give the most consistant results. In all cases the final selection 
for the time for end of thaw was made after considering the base 
temperature, pore pressure and deflection curve. The difficulty 
in pinpointing the time at which the thaw zone reaches the base of 
the ee is felt to be the major cause of the scatter of points 

is 


Smax 


plane reaches the base could be more accurately determined if two or 


on tne 


versus log R diagram. It is felt that the time the thaw 


three temperature sensors were placed at the base of the soil. This 
would provide a more definite indication of the exact time the thaw 


zone reaches the base and should improve the results obtained. 


It will be noted that at the start of thawing, the slope of 
many of the settlement curves increases until a constant value is 
obtained after five or ten minutes. This is because several minutes 
are required to come up to the desired temperature on the top load cap. 
Until a constant temperature is reached the rate of thaw is always 
increasing. Furthermore, there is some water in the top porous stone 
which must melt before the thaw plane proceeds into the soil and this 


thaw would be accompanied by very little settlement. 


It is possible to compare the experimental values for a with 
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values computed from the thegretical Neumann solution 
where the value for a is given by: 
a = 2) Veo (5.6) 


where Ko denoted the diffusivity of the thawed region, ~ 


a 


and -» denotes the root of the following equation: 

-\2 k K V2 T 4G Ky/Ko xh ee 
aa RQ GD ee? ce 
erfr k K V-T eos, AWA CAV AT. 

2 i. 1 2 2 ] 
erfc [(—) Al 
a 
(527) 


where the subscript '1' refers to the solid phase and the 

Subscript '2' refers to the liquid phase, 

kK denotes the conductivity, 

«k denotes the diffusivity, 

Ty denotes the difference between the initial temperature 
of the solid and the freezing temperature of the liquid, 

VY denotes the difference between the initial temperature of 
the solid and the step temperature applied to the Surface, 

C denotes the volumetric specific heat, 

be denotes the volumetric latent heat of fusion of the soil, 


erf (x) denotes the error function, 


and ertce(x) =9108 erf(x). 


The various thermal properties for the soil can be determined 


from information given by Kersten (1940). 


The specific heat is found from: 


C (0.24% 095%) (5.8) 
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and C, Fae (0.2 + w) (5.9) 
where ee denotes the dry density of the soil, 
and w denotes the moisture content. 


The conductivities Ky and Ko can be found for the sojl 


from charts provided by Kersten. 


The diffusivities Ky and «x, are found from 


2 
peas! 
Late (5.10) 
and 
2 ©, (5.11) 


The volumetric latent heat of fusion for the soil can be 
determined from 


Le = Oe oie} 
where L denotes the latent heat of fusion for water. 


The initial temperature and step temperature are known and 
equation (5.7) can be solved for 4. A value for a can therefore be 


determined from equation (5.6). 


The theoretical values for a were computed in this way for 
all the tests carried out on Athabasca clay. The results are presented 
in Table Conl a The theoretical values for a have been plotted against 
the actual values obtained in Fig. (5.8). The figure demonstrates 
that there is excellent agreement between theoretical and experimental 
values. The scatter of the points could be due to experimental error 


in the determination of a or possibly due to the simplifying assumptions 
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Fig. (5.8) Relationship between Rates of Thaw (a) determined 
Theoretically and Experimentally (Athabasca Clay) 
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that are made in determining the thermal properties for the soil in 


the theoretical solution. 


The thaw consolidation ratio is calculated from the coefficient 
of consolidation which is determined from the post thaw settlement 
cunve. [ft Cy, has a substantially different value during the thaw 
portion of the test as compared to the measured post thaw value, 
then the value obtained for R will be in error. The sername of Cy 
with void ratio for the tests on Athabasca clay is shown on Fig. (5.4). 
It will be observed that for the majority of tests Cy varied from 


Sh pe ; . : ; 
cm/sec while the average void ratio ranged in 


Oem tons. 02 cn Om 
value from 0.65 to 1.1. Fig. (5.4) shows that the greatest variation 
in Cy occur at lower void ratios. The pore pressures obtained from 
these tests indicate that R was correctly calculated using the C\, 
obtained from post thaw settlement data. nanan tests should 
be carried out to ascertain the benavior of Cy at very high void 
ratios. Indications are, however, that provided a very large change 
in void ratio does not occur during the thaw phase of a test, then 


the Cy obtained from the post thaw settlement curve will provide an 


essentially correct value for the overall test. 


There is some error introduced when the final height of 
the sample is measured after it is removed from the permode. This 
error will affect the thaw-consolidation calculations in a negligible 
way and will cause the e log P curve to be shifted slightly up or 
down. The quality of the sample obtained by using the method of 
preparation previously described was found to be excellent. Visual 


examination showed the samples to be completely homogeneous throughout 
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and to contain no visible air voids. The final samples were found 
to be 100% saturated with the exception of the sample removed from 
the test 5-5. The latter was 94% saturated which may be due to the 
presence of air or to inaccuracy in the final height determination 


for this sample. 
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9.3 Recommendations for Improvements to Apparatus 


The permode can be substantially improved in its design in 
order to facilitate its use and improve results. It is very difficult 
to set up the existing apparatus for tests on undisturbed frozen 
samples. The present pore pressure measuring system is prone to 
failure through leaks or plugging of the tubes. The outer steel 
jacket is not required for the low pressures encountered in testing. 
Provision should be made to measure deflection from a point located 


directly on the load cap. 


One improved design which eliminates most of the above 
problems is shown in Fig. (5.9). The incorporation of a closed back- 
pressure consolidation cell into the permode has not been carried 
out in the new design. It is felt that the advantages to be gained 
through the use of back pressuring are outweighed by the increased 
complexity of the system. Back pressuring will be required in order 
to investigate the pressures generated during the freezing process 
because of the magnitude of the pressures involved. It would be 
preferable to develop a separate specialized piece of equipment 
to perform these investigations once the specific nature of the 


proposed experiments has been ascertained. 


A closed back pressured system would permit the 
application of high loads on the sample through the use of air pressure. 
The loads required for most thaw-consolidation testing are less than 
5.0kg/cm°. This pressure range can be easily attained with a direct 


loading system. The direct loading system also allows for 
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counterbalancing of the load hanger to allow for friction of the 
load cap. Furthermore, the pressure which can be applied to the load 


cap is limited by the compressive strength of the thermoelectric 


element (10kg/cm? on the sample). 


A back pressured system permits the sample to be saturated 
to 100% by driving the air into solution under high pressure. It 
has been demonstrated that remolded samples can be placed in the 


permode to achieve very nearly 100% saturation without back pressure. 


Undisturbed samples should be back pressured to simulate 
the pressures which existed in situe prior to testing. This is a 
great advantage of back pressuring but against this advantage must 
be weighed the disadvantage that it is no longer possible to quickly 
and conveniently freeze the sample from the top down by providing 


a heat sink such as dry ice above the load cap. 


Improvements can also be made to the temperature controllers. 
A larger or better quality transistor should be provided to replace 
the present 2N5437 as this transistor is just barely adequate to 
handle the loads. The Variac permits the variation of the input current 
and should be retained. The feedback system should be modified so 
that it will automatically vary the polarity of the D.C. output as 
well as the current. Under the present system if a temperature of 
-10°C is selected and the heat sink is colder than -10°C, the controller 
will prevent the temperatures from rising above -10°C but will not 


prevent them from dropping below -10°C. 
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Pete CONG US LENS 


The data clearly demonstrates that within the framework 
of their assumptions, the physics of the thaw-consolidation problem 
have been correctly evaluated in the theory put forward by 


Morgenstern and Nixon (1971). 


The permode proved capable of simulating the required 
temperature and stress boundary conditions for one-dimensional 


thaw consolidation. 


The validity of applying the Neumann solution to predict 
the movement of the thaw plane for the one-dimensional case has 


been confirmed. 


The fundamental importance of the thaw consolidation ratio 
as a parameter indicating the pore pressures developed at the thaw 
plane has been demonstrated. The rate of settlement during thaw and 
after thawing has stopped has also been predicted correctly by the 
theory. The degree of consolidation which occurs during thawing 
was shown to be dependent on R, however, there was some scatter of 
the experimental points. This is primarily due to inaccuracies 
in measuring deflection and the difficulty in determining the exact 


time that the thaw plane reached the base of the sample. 


The thaw consolidation ratio R will provide a useful field 
guide in determining which soils will become unstable under given 
applied pressures and thermal loadings. Field evidence (Brown and 


Johnson, 1970) suggest that in many instances thawing does occur almost 
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proportional to the square root of time. Where thawing does not occur 
proportional to the square root of time, it can be shown that the 
parameter R is not independent of time. Investigations should therefore 
be carried out to determine whether or not the assumption of thawing 
proportional to the square root of time will cause significant 
differences in predicted pore pressures and settlements when other 


rates of thaw occur. 


The pore pressures obtained were generally slightly higher 
than the pressures predicted by the theory. Indications are that 
this effect is due to non-linearity of the e - o' relationship. If 
this is the case, then the effect would be more pronounced for soils 
at hign initial void ratios. Further investigation is required 


in this area. 


The assumption that Cys for the overall test, can be 
determined from the post tnaw settlement curve appears to be acceptable 
for the range of void ratios tested. Investigations should be carried 
out to determine the behavior of Cy at high void ratios and hence 
determine the validity of this assumption for soils at high initial 


void ratios. 


The. initial effective stress one) was assumed to be zero 
for the tests conducted and the results indicate that this assumption 
was essentially valid. It is recognized, however, that under certain 
circumstances tne initial effective stress could attain a Significantly 
high value and the pore pressures obtained would be much less than 
the predicted values. An investigation into the degree of overconsolid- 


ation which occurs as a result of the freezing process should be made 
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and the results obtained related to the pressures which occur on thaw. 
Such an investigation would serve to further the present understanding 


of the whole freeze-thaw phenomenon as related to porous media. 
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APPENDIX A 


SUMMARY OF SOIL PROPERTIES 
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ATHABASCA CLAY 


Unified Soil Classification CL 


Liquid Limit 40.8 
Plastic Limit 2083 
Plasticity Index 20.5 


G, 205 


Approximate Consolidation Parameters (No Freeze Cycle) 


4 


Gp > Mee 1097 “to Mae x 1073 an-/sec 


k= ex 10°? to. 18x Tame cm/sec 


Grain Size Curve Attached, See Fig. (A.1). 
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LAKE EDMONTON CLAY™ 


Unified Soil Classification CH 


Liquid Limit 70.0 
Plastic ime 2 Av 
Plasticity indexs41).0 


G. 2./6 


Approximate Consolidation Parameters (No Freeze Cycle) 


os 1X 1074 one ece 


KR = Ay to 1072 cm/sec 


Grain Size Curve Attached, See Fig. (A.2) 


* 
Source 

Thomson, S., 1969, "A Summary of Laboratory Results on 
Lake Edmonton Clay" Unpublished Internal Note: SM2, University of 


Alberta. 
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BENTONITE” 


Unified Soil Classification CH 


Liquid Limit 59.11. 
Pl dS Cicada mat 87d. 
Plasticity Index 504. 


Bhs Za 


Approximate Consolidation Parameters (No Freeze Cycle) 


eae Ourax 107° ena ieee 


eo ae eee 1072 cm/sec 


Grain Size Curve Attached, See Fig. (A.3). 


* 
Source 


Ganapathy, G.V., 1969, "Some Fundamental Aspects of the 


Residual Strength of Clays", Unpublished M. Eng. Report, University 
of Alberta. 
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APPENDIX B 


SPECIFICATIONS OF MAJOR COMPONENTS 
OF 
TEMPERATURE CONTROL DEVICES 
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THERMOELECTRIC ELEMENTS 


MC5-8-12 Pre-insulated surface, Thermo-electric modules. 
Available from Materials Electronic Products Corporation 


990 Spruce Street, Trenton, N.J. 08638 


THERMISTORS 


YST 511 Non-interchangeable thermistor probes. 


Available from Yellow Springs Instrument Co., Incorporated, 


Yellow Springs, Ohio. 


TEMPERATURE RECORDER 


Honeywell Electronik 15 Strip Chart Multipoint 
Recorder 15303836-24-04-0-000-004 10 
Available from Honeywell Industrial Products Group, 


Wayne and Windrim Avenue, Philadelphia 44, Pa. 


TEMPERATURE CONTROLLER - MAJOR COMPONENTS 


VARTAC 
Type W5MT3A Autotrans former 


Available from General Radio Company Concord,Mass. 


Hammond 60HZ Filament Transformers, Model 167V12 
Potentiometers, Model TRC 8400, 20KOHMS 
RD 412 IRC Dials 

Available from Cardinal Industrial Electronics Ltd., 


14310 - 118th Avenue, Edmonton, Alberta 
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Transistors, Germanium Power, PNP, Type Mota ZN5437 
Transistors, Silicon Power, NPN, Type Mota ZN5067 
Transistors, Silicon Power, PNP, Type Mota ZN56/9 


Available from Bowtek Electric, Vancouver, B.C. 
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SUMMARY OF DATA 
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APPENDIX E 
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Test Series 5 - 6 
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